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Abstract 

We present the Mathematica package HypExp which allows to expand 
hypergeometric functions jFj-i around integer parameters to arbitrary 
order. At this, we apply two methods, the first one being based on an 
integral representation, the second one on the nested sums approach. The 
expansion works for both symbolic argument z and unit argument. We 
also implemented new classes of integrals that appear in the first method 
and that are, in part, yet unknown to Mathematica. 



Package summary 

Title of the Package HypExp 
Version 1.0 

Package obtained from 

http : //www-theorie . physik . unizh . ch/~inaitreda/HypExp/ 

E-Mail maitreda@physik.unizh.cli, thuber@physik.unizh.ch 

Licence none 

Computers Computers running Mathematica under Linux or Windows 
Operating system Linux, Windows 
Program language Mathematica 

Memory required to execute depending on the complexity of the problem 
Other Package needed the package HPL, included in the distribution 
Externed file required none. 
Keywords Hypergeometric functions, expansion 

Nature of the physical problem Expansion of hypergeometric functions 

around integer- valued parameters. These are needed in the context of di- 
mensional regularisation for loop and phase space integrals. 

Method of solution Algebraic manipulation of nested sums and integral rep- 
resentation. 

Restrictions on complexity of the problem Limited by the memory avail- 
able 

Typical running time Strongly depending on the problem and the availabil- 
ity of libraries. 
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1 Introduction 



As solutions of a large class of differential equations, hypergeometric functions 
pFq appear in many branches of science. They appear, in particular, in par- 
ticle physics during the calculation of radiative corrections to scattering cross 
sections in loop [1-6] or phase space [7-9] integrals. In the context of dimen- 
sional regularisation, the arguments of the hypergeometric functions have to be 
expanded in a small parameter around integer or rational arguments. 

Until recently, there was no systematic approach to the expansion of hy- 
pergeometric functions. The required expansions have been produced with a 
case-by-case approach. Recently a general algorithm has been developed [10] 
for expanding hypergeometric functions and other transcendental functions sys- 
tematically around their parameters. This algorithm was implemented [11] in 
the framework of GiNaC [12]. Related work was also presented in Refs. [13,14]. 

Until now, an implementation of the expansion of hypergeometric functions 
around their parameters was missing in the widely used computer algebra sys- 
tems Mathematica [15] and Maple [16], allowing to use hypergeometric functions 
in connection with the multi-purpose features of these programs. 

The aim of this work is to provide an implementation of these expansions for 
Mathematica. This implementation is prepared in the form of a Mathematica 
package that aims to be tunable and user-friendly. 

This paper is structured as follows. In the next two sections, we present the 
two approaches we used in the package. Sectionals devoted to hypergeometric 
functions of unit argument. In section |S1 we describe the implementation of our 
methods in the package HypExp and provide examples of its usage. Finally, we 
conclude with a summary. 



2 Method of integration 

Among the hypergeometric functions (HF), the ordinary GauB-hypergeometric 
function 

2Fi(Ai,A2;Bi;z) (1) 

appears most frequently in scientific calculations, and it is therefore worth to fo- 
cus in particular on this type of functions [17]. We will calculate the e-expansion 
of the 2i^i-functions up to and including order 0(6^) by means of the algorithm 
described below which is based on the well-known integral representation of 
the 2^1-functions. Applying this method to this subset of functions has by all 
means its benefits since it is faster and more efficient than the later described 
nested sums method. However, going to higher orders in the e-expansion or to 
higher j_Fj_i-functions quickly reveals that the nested sums method is in gen- 
eral more powerful. Nevertheless, applying several independent methods also 
provides useful consistency checks. 

2.1 2^1-algorithm 

We first set up our notation. In this work we consider 

jFj_, ({Ai, . . . , Aj}; {B,, . . . , z) (2) 

with 

Ai ~ ai + ait, Bi = bi + /J^e, a^, bi £ Wj and a^, Pi G M. (3) 

We start the description of our algorithm by defining some subsets of C. Let 

V := {zeR\z>l} (4) 
W := €\V. (5) 
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Furthermore, we state that the mapping 

/W:=J^ (6) 
as a Moebius transformation is a bijective mapping and satisfies 

f{V) = V f{W) = W, and f{f{z))=z. (7) 

We finally turn our attention to the 2-^1 -functions and will first collect some 
useful formulas [17, 18]. We note that we can shift each of the parameters up or 
down by integer units via 

,F^{A,,A,-Br,z) = ' + ' " .F,{A, + l,Ar,B,-z) 

Ai — til + I 

A2 — Hi + I 

, + (9) 



{A2 ~Bi + l) 
B1-A1-A2- 

Bi{z-: 

{Bi-Ai + l)[Bi-A2 + l)z 



2Fi{Ai,A2;B,;z) = ^lEl^A^-^l±ni^ ,Fi{Ai, A2;B, + 1; z) 

Bi[z - 1) 



Bi{Bi + l){z - 1) 



,Fi(Ai,A2;Bi + 2;z), (10) 



2FiiA„A2;Bi-z) = " ' " ^^^^^^^ ' ^' ^''^ 

+ 2Fi{Ai,A2 -2-Bi-z), (12) 

- (^i-l)[2-i?i-(^i+^2-2i?i+%] 

2-f'l(,^l, ^2, Hi, — — 



{Ai-Bi + l){A2-Bi + l)z 
X2Fi(Ai,7l2:Bi - l:z) 

{Bi-l)(Bi-2)[z-l) 



{Ai-Bi + l){A2-Bi + l)z 
These relations all stem from the formulas 



2i^i(Ai,A2;Bi -2;z). (13) 



{Bi-Ai)2Fi{Ai~l,A2;Bi;z) 

+ [2A1-B1 - {Ai-A2)z\ 2Fi{Ai,A2;Bi-z) 

+ Ai{z- 1) 2i=^i(^i + 1, A2; Bi-z) = Q, (14) 

{Bi-A2)2Fi{Ai,A2-l;Bi-z) 

+ [2^2 - Si - (^2 -Ai)z] 2^1(^1,^2; Bi; ^) 

+ ^2(^ - 1) 2-Fi(Ai, A2 + 1; Bi-z) = , (15) 

Bi{Bi - l)(z - 1) 2Fi(Ai, A2; Bi - 1; 2) 

+ 5i [Si - 1 - (2Bi - ^1 - A2 - l)z] 2^1(^1, A2; Bi- z) 

+ {Bi - Ai){Bi -A2)z2Fi{Ai,A2 + l;Bi + l;z) = 0. (16) 
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There is yet another class of relations between 2^'i-functions, namely the 
relations of GauB between contiguous functions [17]. Their inclusion would lead 
only to a minor simplicifation here, and thus we let them serve as a check for 
our results rather than implementing them in our algorithm. 

From the relations © - (fT!^ we conclude that the knowledge of the e- 
expansions of the 2-^1 -functions whose integer parts {ai, 02, 61} of the parame- 
ters read 

{0,0,0} {0,1,0} {0,0,1} 

(17) 

{1,1,0} {0,1,1} {1,1,1} 

is sufficient in order to derive the e-expansion of any 2fi-function with {oi, 02, 61} 
being integer- valued. But even this small set of functions can be reduced further 
by means of Kummer relations [17,18]. The relevant Kummer relations read 

2^1(^1,^2; = il-zf'-^'-^-'2FiiBi~Ai,Bi^A2;Bi;z), (18) 

2Fi{AuA2;Bi;z) = (1 - z)-^i 2-^1(^1, Si - A2; Si; ^-) , (19) 

z — 1 

2^1(^1,^2; Bi;z) = (l-z)-^^2i^i(Bi-^i,^2;Si;-^). (20) 

and relate both the functions {0, 1, 1} and {1, 1, 1} to the function {0, 0, 1}, such 
that we can get along with a mere four functions, namely 

{0,0,0} {0,1,0} {0,0,1} {1,1,0}. (21) 

For completeness, we mention that Eq. (|18|l holds true for all z G C, whereas 
Eqs. ifT^ and (EOJ are only valid for z £W. 

The sets (jl7|l and H21|l of basic hypergeometric functions have, however, one 
major drawback. In order to express a general 2-F'i-function solely in terms of 
functions from these sets, repeated application of Eqs. (jHl - (tT^ is required 
and additional negative powers of e might be generated in prefactors via this 
procedure. It is therefore necessary to know the e-expansions of the basic hy- 
pergeometric functions to higher order than is sought by the 2^1-function in 
question. 

In order to avoid this disturbing feature we consider an extended set of basic 
hypergeometric functions. The extended set has three major subsets. In the 
first subset we collect those basic HF's that contain only positive integer parts, 
namely 

{0,0,0} {0,1,0} {0,0,1} {1,1,0} 

(22) 

{0,1,1} {1,1,1} {0,1,2} {1,1,2}. 
The second subset contains those basic HF's in which bi = 0. It reads 

{-1,-1,0} {-1,0,0} {-1,1,0}. (23) 

The third subset finally contains those basic HF's in which bi = —1: 

{-2,-2,-1} {-2,-1,-1} {-2,0,-1} 

{-2,1,-1} {-1,-1,-1} {-1,0,-1} (24) 

{-1,1,-1} {0,0,-1} {0,1,-1} {1,1,-1}. 

Although some functions in this set might not be considered independent since 
they are related via Kummer relations, we will consider this set as basic since 
it will allow us to implement the algorithm described below efficiently and con- 
veniently. The goal of the latter is to express a general 2-F'i-function entirely in 
terms of functions from the set (|22|l - l|24(l by repeated application of both the 
equations (jHl - l|13|l and the Kummer relations (|18|l - (|2()|l . Before we start, we 
mention that throughout the algorithm we make use of the symmetry Ai ^ A2 
after each step in order to ensure that we always have ai < 02. 
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1. We start the reduction of our 2i^i-function in question by applying Kum- 
mer relations such that the sum 



gets minimized. Especially for high absolute values of the parameters this 
procedure shortens the algorithm significantly. 

2. Then, if bi < -1, we apply ifTIHl repeatedly to all HF's with bi < -1. 
This step ensures that from now on we only have to deal with functions 
in which bi > —1. 

For the rest of the algorithm wc distinguish two cases, namely hi — —1 and 
bi > 0. The further steps for 6i = — 1 and 6i > are illustrated by the flow- 
charts in Figures n and El respectively. The ambitious reader is invited to verify 
that at the end of this algorithm only 2^i-functions from the set H22(l - H24|l 
appear and that no negative power of e has been generated at any intermediate 
step. To conclude this section we remark that for the most frequent case in 
which the three parameters oi, 02, and bi are all non- negative, only the last 
column of Figure [3 has to be considered. 

2.2 Expansion of the basic 2^1-functions 

Now that we went through the algorithm for the 2-P'i-functions in great detail 
we have to explain how the e-expansions of the basic HF's from the set H22|l - 
(|24|l are obtained. 

For the functions in Eq. (|22|) we adopt the integral representation [17] 



which we must restrict to z ^ W and Bi > A2 > 0. Since this approach will 
be based on the requirement that the integration over u and the expansion 
in e commute, we have to set up the additional condition &i > a2 > 0. One 
recognizes immediately that from the set H22f) only the functions {0,1,2} and 
{1,1,2} satisfy the latter inequality. The parameters of the other six functions 
first have to be shifted by means of Eqs. lO and Hl()(l until a convergent integral 
representation is obtained for each of them. The subsequent expansion of the 
integral representation in e and how one solves the occurring integrals is covered 
in section 

The functions in Eqs. (|23() and (|24|l . all of which contain at least one nega- 
tive parameter, are now, for the sake of obtaining their e-expansion, expressed 
in terms of functions from the set 122|l . This is again done by appropriate 
application of the Kummer relations (|18|l - H2()(l as well as Eqs. (|Hl - 110(1 . 

The shift of parameters as described in the preceding two paragraphs is 
now unavoidably accompanied by the advent of negative powers of e in certain 
prefactors. To be more precise, we must expand the functions {1,1,0} and 
{1, 1, 2} up to and including order ©(e^) and the other six functions of Eq. H22|l 
to order 0(6^). This might at first glance seem peculiar, but it turns out that at 
the respective highest order in e integrals of the same type appear. We conclude 
from this that, with the tools provided here, it is in principle possible to expand 
a certain class of 2-F'i-functions even to order 0(e^). However, for reasons of 
simplicity and clearness we apply the method of integration throughout only for 
the expansion of 2-Fi-functions up to and including order 0{e^). 



\ai\ + \a2\ + |6i| 



(25) 




(26) 
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a2 = -2 

and 
ai < -2 



apply © 
repeatedly to 
all HF's with 

ai < -2 






{-2,- 
{-1,- 


2,-1} 
2,-1} 



bi = -l 



a2 < 
ar 

ai < 


; -2 
id 

; -2 






apply © 
repeatedly to 
all HF's with 

a2 < -2 



^ {ai,-l,-l} 



L {ai,-2,-l} 



«2 > -1 

and 
ai < -1 



apply dHl) 
repeatedly to 
all HF's with 

ai < -2 



{- 


-2 


-1, 


-1} 


{- 


-1 


-1, 


-1} 




{- 


-2, 


>o, 


-1} 


{- 


-1, 


>o, 


-1} 



a2 > -1 

and 
ai = —1 



apply m 
repeatedly to 
all HF's with 

(32 > 1 



{-1,-1,-1} 

{-1,0,-1} 

{-1,1,-1} 



apply flT^ repeatedly to all HF's with a2 > I 



{-2,0,-1}, {-2,1,-1}, {-1,0,-1}, {-1,1,-1} 




apply (HH) 
repeatedly to 
all HF's with 
ai > 1 



apply (HI 
repeatedly to 
all HF's with 
a2 > 1 



{0,0, 


-1} 


{0,1, 


-1} 


{1,0, 


-1} 


{1,1, 


-1} 



Figure 1: 2-F'i-algorithm for bi = — 1 



02 < —1 and tti < — 1 



apply repeatedly to all HF's with 02 < —1 












i 


{ai,-l,0} 




{ai,-l,> 0} 


{ai,0,> 0} 











apply (0) apply © 

repeatedly to repeatedly to 

all HF's with all HF's with 

Oi < — 1 ai < —1 



{-1,-1,0} 
{0,-1,0} 



-\ {0,-1, >0} - 
^{-1,-1,>0} 



apply (jn} once 



{-1,0,>0}, {-1,1,>0} 



61 > 



02 > and ai < 



apply dS} repeatedly to all HF's with ai < — 1 



£ 



1 




{-1,>0,0} {-1,>0,>0} {0,>0,>0} 



apply (Hg) 
repeatedly to 
all HF's with 

(32 > 1 



{-1,1,0} 
{-1,0,0} 



apply I 
once 



{0,>0,>0} 
{1,>0,>0} 



apply (HH) 
repeatedly to 
all HF's with 
ai > 1 



apply (HI 
repeatedly to 
all HF's with 
02 > 1 



{0,0,>0} 
{0,1,>0} 
{1,0,>0} 
{1,1,>0} 

1 



apply ()13|) repeatedly to all HF's with 
both (61 > ai + 1) and {bi > 02 + 1) 

1 



functions from the set ()22|1 



ure 2: 2 -Fi- algorithm for 61 > 



2.3 Integrals and their algorithm 



2.3.1 General aspects 

Expressing a hypergeometric function 2^1 (ii + cuie, a2 + Q;2e; hi + (3it] z) with 
&i > a2 > in terms of its integral representation according to Eq. H26|l and 
subsequently expanding in e yields integrals that have the general form 

/(X1,X2,X3,X4,X5,^) := /dw , , ^ (27) 



where the Xi sre non-negative integers and z £ W . The results of these inte- 
grals contain rational functions, logarithms, poly logarithms Lz„, Nielsen poly- 
logarithms Sn,p, and Harmonic polylogarithms -ffmi,....mfc, a-U of which will be 
explained in more detail in sections 13 . 1 . II and in appendix lA.il 
Let us define the weight w of the integral by 

X2 + X3 +X4 + 1 - (5x5,0 ■ (28) 

The weight w is related to the weight of the harmonic polylogarithms defined 
later, namely any integral of weight w can be expressed in terms of harmonic 
polylogarithms of weight w or less. In order to guarantee that the e-expansion 
of any 2^'i-function up to order ©(e") can be performed, the computation of 
all integrals with weight w up to rt + 1 is required. 



2.3.2 Description of the algorithm 

In order to make the computation more efficient we now show that any integral 
I (Xij X2, X3j X4: X5j ^) can be expressed in terms of integrals with Xi = X5 = 
and of at most the same weight as the original one: 

1. In the case xi < Xs, the integral / (xi, X2, Xa, X4, X5, 2;) can, via the re- 
placement 

u^-{uz-l) + - (29) 
z z 

in the numerator and subsequent cancellation of the corresponding de- 
nominators, be expressed as a linear combination of integrals of the form 

I (0, X2, Xa, X4, X5,z), I (0, X2, Xa, X4, Xs 

■■■ , I{0,X2,X3,X4,X5-Xi,z) ■ (30) 

2. In the case xi > XSi we can apply the same steps as before and eventually 
arrive at integrals of the form 

I (Xi - X5, X2, Xa, X4, 0,2;),/ (xi - X5 - 1, X2, Xa, X4, 0, z) , 
■■■ , I{0, X2, Xa, X4, 0,2;),/ (0, X2, Xa, X4, 1, z) , 

••• , ^(0,X2,Xa,X4,X5,2;). (31) 

In other words, the replacement rule (|29|l for u yields integrals in which 
at least one of the entries xi and xs is zero. 

3. We continue our reduction by manipulating integrals of the form 
I (Xi, X2, Xa, X4, 0, z) with xi > 0. By recalling the relation 

1 

fdu4- (1 - u) In'^^u) In^^l - u) ln^\l - zu)] = (32) 
J du 
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and explicitly taking the derivative of the product, one derives the relation 

(Xi + 1) ^ {Xi,X2, X3, X4, 0, z) = 

Xi I (Xi - 1, X2, X3, X4, 0, z)~X2l (Xl, X2 - 1, X3, X4, 0, z) 
+X2 I (Xi - 1, X2 - 1, X3, X4, 0, z) - X3 (Xl, X2, X3 - 1, X4, 0, z) 
+ZX'iI {Xl,X2,X3,X4 - 1, ^,z) - ZXiI{Xl + 1,X2,X3,X4 - 1,1, z) . 

(33) 

Repeated application of steps^ — 121 finally yields an expression which contains 
only integrals with Xi = ^'^d of at most the same weight as the integral we 
started with. 

4. The remaining task is now to subsequently lower xs- By repeated appli- 
cation of the relations 

X5-1 



^(0,X2,X3,X4,X5,2;) - X! \ \k] ^('^'^2,X3,X4,fc,z) 

fc=0 ^ ^ 

(X5-1)! h^ixA-v + i)\[j^^ \ J j"- ^ ^ \ 

X ^/(0,X2,X3,X4-t^ + l,0,z)} (34) 



for Xs > 1 and 

I (0, X2, X3, X4, 1, z) = -I (0, X2, X3, X4, 0, z) 

z d 



/(0,X2,X3,X4 + l,0,z) (35) 



X4 + 1 dz 

one eventually arrives at an expression of the desired form. 

5. We can, however, reduce the number of distinct integrals even further. The 
transformation u — > 1 — w in the integrand of H27I) allows us to represent 
integrals with X2 < X3 in terms of integrals with X2 > X3- -By means of 
the relation 

In [1 - z(l - u)] = In (1 - z) + ln(l — u) , (36) 

z — 1 

which holds true for all z £ W and < m < 1, one easily derives the 
formula 

I (0, X2, X3, X4, 0, z) = 5] r J ) • In^^-" (1 - z) ■ /(O, X3, X2, m, 0, ^). 

(37) 

Since this formula also transforms the argument z, we derived relations 
between polylogarithms that allow to simplify again these arguments. 
The argument transformations of the polylogarithms are described in sec- 
tion and in appendix I A. II 

To summarize, for a given weight w the set of basic integrals consists of all 

-'"(0,X2,X3,X4,0,z) with X2 + X3 + X4 = u'- (38) 

Integrals thereof with X2 < X3 can be reexpressed in terms of integrals with 
X2 > X3 a-iid argument z/(z — 1) and are therefore easily obtained once the 
argument transformations of the polylogarithms are taken into account. 
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We implemented this algorithm in the function HypExpInt. This function 
will be explained in detail in section 15.11 For the actual calculation of the e- 
expansion of 2-^1 -functions we need a bit less than is provided by the algorithm 
and by the function HypExpInt. We have: 

XI < 2, (39) 
X5 < 1. (40) 

The first inequality arises from the fact that we have to shift the parameters of 
our basic hypergeometric functions, Eq. (|22|l . via the relations and pOfl in 
order to obtain convergent integral representations. In this procedure, bi — 2, 
which eventually determines the highest value for xi to occur, assumes values 
up to 2. Similarly, the parameter oi determines the highest value for ^5 that 
can show up. With this in mind one derives the second inequality directly from 
the collection (|^ of basic hypergeometric functions. 



2.3.3 Integrals of unit argument 

Putting 2 = 1 in the integrals / (xij X2, X3, X4, X5, ^) immediately turns our 
attention to an other type of integrals that we considered useful to implement. 
We define the function U {n,m,p) by 



U{n,m,p):= /du ln"(u) • ln™(l - u) • (41) 





with p £ I1 and n, m being non-negative integers. In the case p < the 
inequality to -I- p > has to be satisfied in order to yield a convergent integral. 
We demonstrate below how any convergent integral U (n, m,p) can be expressed 
in terms of U (0, 0, 0) and integrals of the form U {n, m, —1). 

1. We start by considering the case p < — 1; rt > and m+p > 0. Repeated 
application of 

n 171 

U (n, m, n) = ■ U (n — l,m,p) -\ ■ U (m — 1, n, —1) 

^ p+1 ^ ^ p+1 ^ ' ' ^ 

+ ^ E U{n,m-1,-T) (42) 
p+1 

leaves us with integrals of the form J7 (n, m, — 1) and U (0, to, p), where the 
latter type still happens to have p < — 1 and to + p > 0. These integrals 
get reduced via the recursion relation 

U (0, m,p) = ^ C/ (0, m - 1, -k) (43) 

to integrals with p = — 1 . 

2. We now proceed with the case p = 0, in which the function U is symmetric 
in n <-> TO and therefore n > m can always be achieved. Applying this in 
turn with the formula 

[/(n,TO,0) = -nJ2 I ^' U {n-l,a,0) 

-fn^^^- :;7(n-l,A,-l) (44) 

for {n, to} 7^ {0, 0} eventually yields integrals of the desired form. 
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3. In the case p > we apply the same trick as in step El of section 1^.3. 21 and 
derive from 

1 

fdu— K(l - u) ln"(w) ln™(l - u)] = (45) 
J du 



the recurrence relation 

{p + 1) U {n,m,p) = p J7 (n, TO,p — 1) + ri J7 (?i — 1, m,p — 1) 

—n U {n — l,m,p) ^ mU {n,m — l,p) . (46) 



Repeated application of these steps finally yields an expression which only con- 
tains integrals of the demanded form. 

U(0, 0, 0) is trivially found to be unity. A nice algorithm for the computation 
of integrals of the form U{n, m, —1) is given in section 7.9.5 of Ref. [19] and will 
not be repeated here. 

The evaluation of integrals [/(n, m,p) with p (z and n, m being non- 
negative integers can be called with the function HypExpU, and an example can 
be found in section 1531 

The connection to / {xi,X2, X3, X4, X5,z) in z = 1 is given by 

/(X1,X2,X3,X4,X5,1) = X](-1)'"''' ( ^/ )c^(X3 + X4,X2,J-X5) (47) 

for X2 > X5- -By means of this relation the function HypExpInt can be directly 
called with unit argument, see section [5m 



3 Nested sums method 

3.1 Definitions and auxilary functions 

In this section, we briefly review the S and Z sums and their properties intro- 
duced in [10]. Their definitions are 

Z(n,{mi,...,mk},{xi,...,Xk}) = XI XI ' ' ' 51 ■ • • ' (^8) 

il = ll2 = l ■ifc = l " 

S{n,{mi,...,mk],{xi,...,Xk}) = ' ' ■ ■ ■ ■• (49) 

21—122 — 1 2fc— 1 ^ ^ 

or the equivalent recursive definitions 

" x^^ 

Z(n, {mi,...,mfc},{a;i,...,Xfc}) = ttht^I^i ^ 1, »2,...,fc, a:2,....fc) (50) 



41 = 1 1 



n ii 



S{n,{mi,...,m,k},{xi,...,Xk}) = r^S{ii,i2.,...,k,X2,...,k) (51) 



ii=i 1 

Zin,{}A}) - [I III (52) 
5(n, {},{}) = W ^<o , (53) 



where we introduced the short-hand notation 

TO2,...,fe = {m2, ■ • ■ 
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The number k is called the depth and the sum of the |mi|'s is called the weight 
of the nested sum. The S and Z sums are related, since their expressions 
differ only by the upper summation limits in the recursion relation. Using 
J2i<j — Si<j ! '^^^ convert S into Z sums and vice versa by means of 
recursive application of the identities 

S{n,mi^...,xi^...) = ^ ^ -^5'(i2,TO3,...,a;3,...) 

ll = l »2 = 1 

+S{n,{mi + 1712,1713^,,}, {xiX2,X3^,,,}), (54) 
Z(n;mi,...,a;i,...) = ^ ^ -^^(«2 - 1, "i3,..., 2:3,...) 

41=1 12 = 1 

-Z(n, {mi +TO2,'Ti3,...},{a;iX2,X3,...}), (55) 

and reconstructing Z sums from the a;*/*™ with help of the definition H48|l . The 
Z and S sums form two algebrae, that means that a product of Z sums is 
expressible as a sum of Z sums. Using 

Z{n;mi^,,,^k,xi^,,,^k) x 2'(n, m'^^ ^;,a;'i^^^ _;) 

^^(n - i,m2....,k,X2....M)Z{ii - l,m[ j^,x[ j) 



.^1*1 



'»2 

X 



X! ";^^(*2 ~ l,"ii,...,fc,a:^i,...,fc)^(i2 - l,"i2,..._i,a;2,...,i) 



i2 = l '2 



(x x' ) 



z=l 



recursively and reconstructing Z sums from the x'/z™ with l|48|) proves the 
claim. Likewise, a product of 5 sums is a linear combination of S sums. Since 
products are defined for equal upper summation limit n, it is useful to have 
a relation between sums Z{n, . . . ) of different n's. The recursive use of the 
following formulae allows to change the upper summation boundary of a Z or 5 
sum, and, doing so, to bring all the nested sums to the same upper summation 
limit. This is called " syncronizing the sums" in [10] 

Z{n + c- l,TOi,...,fe,a;i,...,fe) = Z (n - l,mi„„^k,xi,„,^k) 



. — 

3=0 ^ •' ' 



Z{n - 1 + j, TO2,...,fc, a;2,...,fc), (57) 



S[n + c,TOi,...,fe,a;i,...,fc) = S'(n, TOi,...,fe, xi^...^^) 

+ H + ^■)mi ■^('^ + "^2,....fc, X2,...,fc). (58) 

Since the definition of S and Z sums is for denominators of the form i^™, it 
will also be useful to convert sums of the form 

n ^ 

i—i ^ ' 

to the form of the definition, that is, with c = 0. If the depth of the sum is zero, 
we use 

T- ^ = -> T- + 7 ^> 59) 

(i + c)™ X ^ (i + c - 1)™ c™ (n + c)™ 

i — \ 2=1 
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and otherwise we use 



E- r — Z(i~l,mi k,xi k) — 
(i + c)™ ^ '■■■'^ 

1 " 

— > 7 ^ — Z(i — l,mi k,xi k) 



i=l 
n-1 



Z(n ~ l,mi,...,fc,a;i,...,fc). (60) 



(n + c) 



3.1.1 Relations to other functions 



Special cases of Z and S sums are related to other functions. For finite upper 
summation limit we have 



Z(n,mi,...,fe,l,...,l) = Zrm,...,mi,{n) (61) 
S'(n,mi,...,fe, 1, . . . , 1) = Smi,...,Tnkin),mi > (62) 



where the .^mi,...,mfc (") are the Euler-Zagier sums [20,21] and the 5'mi,...,mfc (^^) 
are the harmonic sums [22]. For infinite upper summation limit, we have the 
following identities 



Z(oO,TOi,...^fc,Xl,...,fc) = Lim^,...,,ni{xk, ■ ■ ■ ,Xi) (63) 

Z(oo,77ii^...,fc, 1, . . . , 1) = C("^fc, • ■ • ,nT-i) (64) 
2'(oo,mi,...,fc,{a;,l,...,l}) = i?,„i,...,,„j^ (a:) (65) 
Z(oo, {n + 1. ■ - ^^ , 1 }. {a:, 1, , 1 }) = S'„,p(a;) (66) 
p—i p~i 

Limk....,mi{xk, ■ ■ ■ are the Goncharov multiple polylogarithms [23]. A spe- 
cial case of the Goncharov multiple polylogarithms are the harmonic polyloga- 
rithms (HPL) Hmi,...,mk(x) of Remiddi and Vermaseren [25]. C,(mk., ■ ■ ■ ,toi) are 
the multiple zeta values [24]. The HPL's can be reduced to classical polyloga- 
rithms {Lin{x)) [19] and Nielsen polylogarithms {Sn,p{n)) [26] up to weight 4 
in our case^. Since for hypergeometric functions we only have one free variable, 
we will only use the last two identities. 

The package HypExp uses the package HPL [27] to deal with the harmonic 
polylogarithms. 

Here one can see the usefulness of the Z sums as a connection between the 
Euler-Zagier sums (which appear in the expansion of the F function, as will be 
shown in section r^.l.2|l and the harmonic polylogarithms which will eventually 
appear in the coefficients of the expansion. 



^ Since for the expansion around integer- valued parameters of the HF, there appear only 
HPL's with positive indices. 
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3.1.2 Expansion of the T function 

r functions can be expanded around integer values as follows [10] 

^ oo 



T{a + ae) = T{1 + ae)T{a) 



1 



l + ^^(ae)^-Zl_l(a-l) I (67) 



r(a + ae) T{1 + ae)T{a) 



1 



l + ^(aeyZi^...,l(a-l) 



/ 



T{l + ae)r{a) 



l + ^(_«e)^5i_i(a-l) I (68) 



for a integer and a > 0. For negative (or vanishing) a, one has to use the identity 

xr{x) = T{x + 1) 



so that 



r{-n + ae) 



r(i 



r(l+ae) 



ae 



i-iy 



n 

r(l - ae) 



ae -^-^ —] + ae 



r(n + 1 — ae) 
r(l + ae) (-1)" ^ 



ae r(n+l) 



l + ^(ae)^-5i,...,l(n) 
j=i ^ — ,^ — ■ 



and 



ae 



r(-n + ae)-r(l + ae)n(-^' + "^) 

r(n + 1 - ae) 



ae 



r(l+ae)' ' r(l-ae) 



(69) 



ae 



T{l + ae) 



(_l)nr(n + 1) I 1 + J2i-aeyZ^_^^_^in) | . (70) 



3.2 Description of the algorithm 

The algorithm described in this section is the adaptation of the algorithm of 
type A of [10] for the special case of hypergeometric functions. 
The definition of the hypergeometric series is given by 



jFj_i({Ai, . . . Aj}, {Si, . . . , Bj_i},x) 

^ {A,)i...{Aj)i x' 
£^(Bi),...(B^_i), i\ 



1 + E 



{Ai),...{Aj), 
^(Bi)i...i?j_i), i\ 

^ ^ r(i?o...r(j3.;_i) g 



r{Ai + i)... r{Aj + i) 



r(Ai) . . . r{Aj) ^ r(Si + i)... r(Sj_i + i)r(z + 1) 



(71) 
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Where {a)i — a{a + 1) ... (a + i — 1) is the Pochhammer symbol. We denote the 
coefficients in front of the sum with E and the sum itself with /. The parameters 
A and B are of the form 

Ai = ai + a,;e Bi — bi + Pie Ui, hi '^li and a^, /S,; £ R . 

In order to give the e-expansion (to order n) of the hypergeometric function, 
one has to expand the product of E and / to order n. Let us have a look at the 
required depth of the e-expansion for these factors. 

• For each negative hj, one gets a factor e^^ in E^ and factors e in /, but 
only for the values 1, . . . , — 6j of i. 

• For vanishing 6j's one gets a factor in E but no factors of e in /. 

• For negative a^'s, one gets factors of e in _B and factors eT^ in / for 
i = 1, . . . , —cij, 

• For vanishing aj's, one gets factors of e in _B but no factors in /. 

This means that for the expansion of the hypergeometric function to order n 
we have to compute the expansion of E to order 

nE=n + #(a, < 0) (72) 

and for I to order 

n/ = n + #(6, <0)-#(a, =0). (73) 
Note that the above argument explains why the expansion of 

2Fi{aie,a2e,b + Pe,x), 6>0 
starts with 1 + O(e^), whereas that of (for example) 

2F1 (1 + aie, 1 + a2e, f}e, x) 

has a term. 

We treat first the prefactor 

^ nB^)...nB,^^) 
T{A,)...V{Aj) 

We expand the V functions with the formulae of section l5.1.2l which leads to 

r(l + /3ie)...r(l+/3j_ie) ^,^ 

E = -pr— ^ 7:7— ^/ (^1, ...,Aj,Bi,..., Bj^i,e) . (75 

r(l + aie) . . .r(l + aje) 

We have to expand the product of Z sums appearing in / to the order 

= + #(& < 0) = n + #(a < 0) + #(6 < 0) (76) 

because of the factors for each negative b. Since the factors r(l + (3ie) and 
r(l + aie)~^ will cancel those of the expansion of the F functions of the sum / 
we will only consider 

^ _ r(l + aie)...r(l + aje) r(Bi) . . . F(Bj_i) 



r(l + Pie) . . . F(l + /3j_ie) T{Ai) . . . T^Aj) 
- f{Ai,...,Aj,Bi,...,Bj^i,e). (77) 

The coefficients of the e-expansion of E are Z or S sums with finite upper 
summation limit. These are rational functions that can be computed easily 
using the recursive definition. 
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We now turn to the sum 



oo 



r{Ai + i)...r{Aj + i) 



T{Bi + i) . . .r{Bj^i + i)T{i + 1) 



(78) 



First one makes use of the identity xT{x) = r{x + 1) to bring all r{Aj + i)'s and 
T{Bj + i)'s to the form r(z + £_e) times some rational functions of i. Then one 
expands the T functions with the formulae of section 1!^. 1.21 As one knows that 
there will appear terms with, at worst, as many factors 1/e as there are strictly 
negative Oj 's, one has to expand the F-functions to the order 



n^j ^ni + #(aj < 0) n + #{bj < 0) - #(aj = 0) + #(0^ < 0) (79) 



in order to be able to get / to the required order nj. The result of this expansion 
is a product of 

• the factors r(l + ^e) from the expansion of the T{i + ^e). We use them to 
cancel those of E. From now on, we denote by / the sum / without these 



• a product of e-expansions with coefficient ^i,...,i(« — 1) or 5*1,....! (i — !)• 
These products are treated as above, first converting the S sums into Z 
sums, then expanding the products of Z sums into a sum of single Z sums. 
In this case one cannot calculate numerical values for the Z sums, due to 
the i in the argument. 

• a rational function R{i, e) of i and e from the factors x in the reduction 
a;F(a;) = r(a; + 1). This is reduced by means of expansion into partial 
fractions to a sum of coefficients times single denominators. 

The next step is to bring the factors l/{i + c)" appearing in R{i, e) to the form 
suitable for applying the definition of the Z sums. The method is presented 
in section IXTl 

At this point, one can perform the last summation over i, as described in 
the following section. The result is a sum of Z sums with upper summation 
limit equal to infinity. These can be converted to more common functions, as 
described in section 13. 1.1 1 

To obtain the full expansion, one has to multiply the expansions of E and / 
and to add unity to the result. 

3.2.1 Last summation 

The algorithm reduces the sum / to a sum of terms of the form 



Here we have to distinguish several cases. The general strategy is to perform the 
summation over i until the denominator is positive for each i, then to simplify 
the remainder to the form of the definition of a Z sum. We list below the 
different kinds of terms that can appear and for each of them the way it is 
processed. 

• Z sums of argument i — 1 times one denominator with negative offset 



factors. 




1 



Z(*-l,...). 



+ ae) 



m 
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Here one has to perform the sum up to the j-th term, according to 

oo ^ 

y] T- — —, — ^^(* - 1> ■ • ■ ) 

j — 1 ^ 

= Yn — Z(fc-1,...) 

k+j 

+ > — Z(fc+j-l,...). (80) 

^ (fc + ae)™ V J ' ^ V ; 

One then expands the first two terms to the required order in e. The 
occurring Z sums can be evaluated, as their upper summation Umits are 
finite. We now take a closer look at the last sum 



oo 



k=l \l=l 



Here one only keeps terms up to the required order in e in the sum over 
The Z sum in the last term has to be syncronized down to argument 
i — 1, as described in section ITTI 

Z sums of argument i — 1 times one denominator without offset 

^Z(z-1,...). 

If n is positive, following the definition of the Z sums, one gets 
^^Z(*-l,m2,...,fc,{l,...,l}) 

i 

= Z{oo,{n,m2,...,mk},{x,l,...,l}) = Hn,m2,....m,,{x). (82) 

In many cases, the harmonic polylogarithm can be expressed in terms of 
more common functions. 

If n is negative, one has to interchange the first two summations 

oo 

^i™x'2'(i - l,m2,...,fe,a;2,...,fc) 

i=l 

oo i—1 i2 

= ^i^x' ^ T7|j2'(i2 - l,nT-3....,fc,a:3,...,fc) 

i=l 12 = 1 2 

oo ^2 ^ 

= E ^^(*2-l,m3,...,fe,a:3,...,fc) E (83) 



12 = 1 j=i2+l 



The last sum 



i=i2 + l 1=12 + 1 



is a polynomial of degree m in 12. The coefficient of each power of 12 is a 
rational function in x. The of this polynomial can be combined with the 
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i^™^ from the definition of the Z-sum. We can also factor out an 

The remainder is then a polynom in i2 with rational coefficient functions 

of X. The result is then 



i'^x'Z m2,...,fe, X2,...,k) 



i=l 



= a; 2^ — Z(i2 - l,m3,...,fe,a;3,...,fc)2_^Cj(a;)i^. (85) 

i2 = l *2 



We are now left over with a sum of terms of the form 

X] ^^^-A ^ (*2 - 1> '^3,...,fe, a;3,...,fe) ■ 



=1 * 



Recursive use of this formula leads cither to Z sums with a denominator 
with positive m or to a Z sum multiplied by a denominator with negative 
exponent n but with zero depth. For the latter we can apply 

Y,xH"^Z{i -!,{},{}) = J^xH"^ = Li_U^)=(x-] (86) 

1=1 i=l ^ ^ 

which is also a rational function in x. 
• Z sums of argument i — 1 without denominators 

a;'Z(j-l,m2,...,fe,{l,...,l}) 
The summation can be performed by exchanging the first two summations 

oc 

Y,x'Z{i-l,m2,...,k,{l,...,l}) 

oo i— 1 ^ 

= H ?HJ^(*2 - 1, m3,...,fe, {1, . . . , 1}) 

i=l i2 = l 2 
oo ^ oo 

= XI ^^(*2-l,m3,...,fe,{l,...,l}) 

82=1 1=12 + 1 

OO ^ . . - 

Yl :^^(*2 - l,m3,...,fc,{l,...,l})- 
Z(oo,rn2,...,fe,{a;, 1,...,!}) 



i2 = l 

X 

1-x 

X 



x 



1-x 

• Single denominator with negative offset 



Hm„...{x) (87) 



Ex^ 
(i — -i 4- 



{i — j + ae)"* 

Here we perform the sum explicitly for the j first terms 



y ^- = 

^ (i - j + ae)™ 

j — 1 j oo 

^ + 7^+^' y TT^^' (88) 
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we can simplify the last summation if we are only interested in the expan- 
sion to the order n: 

J\! — 1 — X \ ^ — / 

i=0 ^ ' k=l 

= E^i^^(-«^)'iwK^)+0(6"+i) (89) 

Note that this works for positive and negative m. 
• Single denominator with positive offset 



5Z a 



X' 



We use the formula 



^ (i + ?■)'" ^ x-J' ^ (« -I- 7)™ ~ a;-?' I «™ ^ 



l(iz,„(x)-^^) . (90) 



4 Hypergeometric functions of unit argument 

The hypergeometric series 

jF,^, ({^1, . . . , AjY {Si, ... , S.-i}; z) = f: ' ' ■^^'\ ^ (91) 



(i?i)„ . . . (i?j_i)„ n! 



converges in z = 1 only if the condition 

.7-1 J 



^ - ^ A > (92) 



t=l r=l 

is satisfied. This implies that, if 

j-i J 

s = ^ 6t - ^ a, > , (93) 

the expansion in e commutes with the procedure of taking the limit z — > 1, and 
the series expansion of the function jFj_i ({Ai, . . . , Aj}; {Bi, . . . , Sj-i}; 1) 
around e = has a well-defined finite radius of convergence. We therefore 
call the case s > non-critical. 

The case s < 0, on the other hand, will be referred to as critical since this 
case requires more care and additional explanation on its treatment. By means 
of an algorithm based on partial fractions it is possible to express a HF jFj^i 
of unit argument and value s in terms of J — 1 hypergeometric functions jFj-i, 
also of unit argument, but of value s 4- 1 or higher. We outline this procedure 
for the case in which no two of the Bi are equal. Always assuming convergence, 
we start with the series expansion (|91ll of the function 

{{Ai -h 1, . . . , A/ + 1}; {Si -M, . . . , + 1}; 1) (94) 
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and multiply and divide therein by the fraction 



(Ai + n) • . . . • (Aj + n) ^gg^ 



(Bi + n) • . . . • + n) 



The inverse of the above expression gets combined appropriately with the T- 
functions of the Pochhammer symbols, whereas the multiplied one gets ex- 
panded into partial fractions, yielding 

The term linear in n gives back the function (|94|l . which then cancels on both 
sides of the equation. The constant term is proportional to the function 

jFj^^ ({^1, . . .,A.,}- {Bi, . . . , i?,7_i}; 1) (97) 

for which we are seeking and for which we can now solve the equation. Each 
term of the last sum in l|9t)|) contains a hypergeometric function of value s + 1. 
After some intermediate steps, we arrive at 



J -I J 

E^.-E^. 



jFj_i({Ai,...,^j};{Bi,...,Bj_i};l) 



i^iI{ilix^r{Br-B^) Br 
XjFj^i {{Ai, . . . , AjY {Si, . . . , Br-uBr + 1, Br+i, . . . , 1) . (98) 

The case in which there is any combination of equal BiS can be treated anal- 
ogously. The linear and the constant term in the expression (|96|l remain un- 
changed, only the last sum will look different, and the final expression (|98|l will 
contain HP's of value s -I- 1 or higher. 

Repeated application of this algorithm allows to express a critical HF as a 
linear combination of non-critical ones, the two expressions being related via 
analytic continuation. 

Since the analytic continuation is unique and the expression obtained by 
the algorithm has a well-defined expansion around e = 0, we can associate the 
expansion at hand also with the original critical HF. In this sense the e-expansion 
of a critical HF has to be understood and the user must be aware of this feature. 
The same phenomenon happens, by the way, for the well-known F-function. 

To elaborate a bit more on this we state that the above procedure works 
only for J > 2 and that the case J = 2 is simpler: 

... X T{Bi)T{Bi- Ai- A2) , , 

-^^(^^'^-^^^^^)- FU-V)r(i?!-Al 

Again, the series expansions of the F-functions have to be understood in the 
sense of analytic continuation. 

The crucial property of a critical function is that one can find at all a value 
for e such that the series (I91II converges, i.e. 

j-i ,7 

^A-^a, ^0 (100) 

in the case s < 0. If s < and Ij 100(1 yields zero, the expression is divergent for 
all values of e and no remedy can be found. The way in which the HFs of unit 
argument are implemented in the package is described in section 15. II 
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5 The Mathematica package HypExp 



We implemented the results of the preceding sections in the Mathematica pack- 
age HypExp. It allows to expand arbitrary -functions to arbitrary order 
in a small quantity around integer parameters and can be obtained from [28]. 
The results are displayed in terms of rational functions, logarithms, polyloga- 
rithms, Nielsen polylogarithms, and harmonic polylogarithms. The results given 
by these functions are not systematically simplified using Simplify, since the 
simplification might take longer than the expansion itself, in particular for ex- 
pansions to high orders. A Simplify might produce a more compact result. 
Furthermore, the package HypExp should be loaded at the beginning of the ses- 
sion. 



5.1 Functions, commands and symbols added 

The package HypExp adds two new symbols 

• $HypExpPath is the path where the HypExp package is installed. 

• $HypExpFailed is the symbol returned by the package in case of failure. 
The package adds the following functions 

• HypExp [Hypergeometric2Fl [ . . . ,x] e,n] gives the e expansion of the en- 
closed hypergeometric function. The function HypExp applied to anything 
else but a HF will leave it intact. Therefore one can map it onto an ex- 
pression containing hypergeometric functions, and only the HF's will be 
expanded to the required order in e. This is illustrated by the following 
example: 



« HypExp ' 

(HypExp[#l, e, 1] S ) 11% 

(Log[l - €] Hypergeometric2Fl [1 + e, 1, 2 

Log L 1 



Log [1 



e] 



6 ( 



Loq[l - x] ^ Log [ 1 



^ PolyLoq[2, X] 



(HypExp[#l, 6, 1] S ) //e 

(Log[l - 6] HypergeometricPFQ[{l + 2 e, 1 
{2 - e, 2 + 3 e), x] ) 



X]) 



2), 



Log[l 



e] (-( 



Log[l 



X] 



6 (-( 



2 Log[l 



X] 



Log[l - x] 
2 X 



2 PolyLog[2, x] 



SiiH)lify [%] 

-((Log[l - 6] (Log[l - X] (2 + 4 6 + e Log[l 
4 E PolyLog[2, X] ) ) / (2 x) ) 



X]) 



The result is not given as a SeriesData since this would have the effect 
of forcing the expansion of the rest of the expression. This example also 
illustrates that the results produced by the package are not simplified, as 
this might be time consuming and not always appropriate. If one wants 
to get a compact result, one should use Simplify or even FullSimplif y. 
The prefactors that accompany the variable e can also be symbolic, and 
the expansion also works for argument 2; = 1 as shown by the following 
example. 



FullSimplify [HypExp [Hypergeonietrlc2Fl [1 + 3 6, 1 
3 + 2 e, 1], e, 3]] 

SeriesData [Global 'e, 0, 



{2, 



-2 (-6 + Pi ) , 36 (-1 + Zeta[3])}, 



HypExp [HypergeometricPFQ[{l + e, 1 - 2 e, 2 
{2 + 2 6, 2 - e), 1] , 6, 2] 

SeriesData [Global 'e, 0, 



, 4, 1] 
3 e). 



<!' I- ^ 



2 Pi 



8 Pi 



78 Zeta[3] 



}, -1, 3, 1] 
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The technicalities of the exansion in the case of the argument being unity 
are explained in sectional 

• The function HypExpInt[xi, X2, X3i X4: X5: A evaluates integrals of the form 

^(Xi,X2,X3,X4,X5,2) with ii; = X2+X3 + X4 + l-'5x5:0 <5, (101) 

as described in section lT^ All the Xi are non-negative integers and z G W. 
The upper bound on the weight w stems from the fact that the expansion 
of 2-^1 -functions up to order C(e'*) also involves integrals of weight 5. The 
integral can be called with the argument being symbolic: 



HypExpInt[2, 0, 1, 1, 


3, z] 


(27 + 2 Pi^ - 24 z - 


4 Pi^ z + 2 Pi^ z^) Log[l - z] 


12 


2 3 
(-1 + z) z 


(-3 + 4 z) Log[l - 


2 3 
z] Log[l - z] 


2 3 
2 (-1 + z) z 


3 z 


2 

Log[l - z] Log[z] 


(-3 + 4 z) PolyLog[2, z] 


2 z 


2 3 
2 (-1 + z) z 


PolyLog[3, 1 - z] 
3 

z 




3 z + 4 Zeta[3] - 


8 z Zeta[3] + 4 z^ Zeta[3] 


4 ( 


2 3 
-1 + z) z 



Arguments that match the pattern z/{z — 1) are treated with the rela- 
tions between polylogarithms of different arguments, section [5.21 and ap- 
pendix being already taken into accoimt: 



HypExpInt [1, 


1- 0. 1, 2, ^ ! ^1 


, (6 + Pi^) 


(-1 + z) ^ Log[l - z] 




2 

6 z 


2 

(-1 + z) 


3 2 2 
Log [ 1 - z] , (-1 + z) Log [ 1 - z] Log [ z ] 


6 


2 ' ,2 
z 2 z 


2 

(-1 + z) 


PolyLog[2, z] 




2 + 

z 


2 

(-1 + z) 


Log[l - z] PolyLog[2, z] 




2 

z 


2 

(-1 + z) 


PolyLog[3, 1 - z] (-1 + z)^ Zeta[3] 


2 2 

z z 



Finally, also all z can be inserted directly. For the special cases z = 
and z = \ the integral simplifies considerably and the restriction w < b 
can be dropped. 



HypExpInt [2, 4, 1, 1, 3, 1] 

3 Pl^ Pi^ 2 

I, + - 42 Zeta[3] + 22 Pi Zeta[3] + 

I Pi^ Zeta[3] - 108 Zeta[3]^ - 348 Zeta[5] + 

16 Pi^ Zeta[5] - 240 Zeta[7] 



In the case z = 1 we refer the reader also to the next paragraph and to 
section l2.3.3l 

• The function HypExpU[n, to,p] is described in section f2.3.3l It evaluates 
integrals of the form 

U {n,m,p) with p G Z and n, m being non-negative integers. (102) 
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The only additional constraint on the parameters is that in the case p < 
the condition m + p > has to be satisfied in order to yield a convergent 
integral. 



HypExpU[4, 3, -2] 






4 Pi^ 
2 PI + ~- - 144 


Zeta[3] + 48 Pi^ 


Zeta[3] + 


18 4 

~ Pi Zeta[3] - 


216 Zeta[3]^ - 


576 Zeta[5] + 


72 Pi^ Zeta[5] - 


1152 Zeta[7] 





• HypExpIsKnownToDrder [ai aj , 61 ,n] returns True if the 

expansion of the hypergeometric function with parameters corresponding 
to the first arguments of the function is available in the library to the 
order n. 

• HypExpAddToLib adds an expansion to the library, this method is de- 
scribed in one of the next sections. 



5.2 Functions modified 

The package also updates Series to allow it to expand compound expressions 
containing hypergeometric functions. The difference between this and the map- 
ping with HypExp is that the other functions of e are also expanded, as shown 
by the following example: 



Series [Log [1 - e] Hypergeometric2Fl [e +1, 1, 2 
{e, 0, 2}] 

SeriesData [Global "e, 0, 



X] , 



, Log[i 



PolyLog [2, x] , 



Log [ 1 - x] 



2 X ' 
1, 3, 1] 



Log[l - x] 



2 



Series [Log [1 - e] HypergeometricPFQ [ { 2 e + 1, 1 
{2-6, 3 6+2), X] , {6, 0, 2)] 

SerlesData [Global "e, 0, 

2 

. Log[l - x] 5 Log [ 1 - x] Log [ 1 - x] 



2), 



2 PolyLog[2, x] 



2 X 

), 1, 3, 1] 



2 X 



This allows to work with the expansion of HF's as Mathematica users are used 
to. We also updated the series expansion of the regularized hypergeometric 
functions since they are nothing else but hypergeometric functions divided by 
F- functions. 

Since the incomplete B function is a special case of HP, 



B{z, a, b) 



■ 2Fi{a,l - b,a+ l,z), a ^ -1, 



(103) 



it is also possible to expand it with the method described in this paper. There- 
fore we also updated the series expansion of the incomplete B function around 
integer values of its parameters, as shown by the following example 



Series [Beta [X, 1+26, 2-6], (6, 0, 1)] 

SerlesData [Global 'e, 0, 



{X - 



2 



-3 X 
2 



2 



^ Log[] 



X Log[l 



4 2 
+ 2 X Log[x] 



- X Log[l 



X Log [x] 1 , 0, 2, 1 ] 



In order to account 



for the relations between polylogarithms Lin and Nielsen polylogarithms Sn^p of 
different arguments we updated PolyLog and implemented the relations (|113|l - 
(|136|l given in appendix I A. II This is also illustrated with an example. 



PolyLog[2, z / (z - 1) ] 

Log[l - z]^ -PolyLog[2, z] 
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5.3 Working with the Ubraries 



Since the computation of the expansion for high orders and large parameters is 

quite time consuming, it is of interest to store the results that have been already 
calculated and reuse them, instead of recalculating them. The HypExp library 
contains expansions of HF for some sets of parameters. When an expansion is 
requested, the package first checks whether the library contains the expansion 
for the requested set of parameters to the requested order. If so, it loads it 
and gives the answer, if not proceeds with the calculation along the line of the 
preceding sections. The library management can be called dynamic in the sense 
that elements of the library are loaded in the memory at run time only if they 
are needed. 

The package HypExp has a standard library. Further libraries can be added, 
depending on the needs and on the amount of available disk space. The different 
libraries can be found at [28]. 

The package also allows to extend the provided libraries with HF's not in- 
cluded in the standard libraries, or included but not to a sufficiently high ex- 
pansion order. The expansion of jFj_i(oi, . . . ,aj,bi, . . . , bj-i,x) to the order 
n can be added to the library with the command 
HypExpAddToLib [oi, . . . ,aj ,bi , . . . ,6j_i,n] 

where oi, . . . , Oj, 61, . . . , are integers. Upon this command, Mathematica 
computes the expansion for arbitrary e-parts added to the integers ai , . . . , Oj 
and &i, . . . , & As arbitrary coefficients are more difficult than numbers, the 
time needed to add an expansion to the library is longer than the time for a 
single evaluation. Therefore adding expansions to the library is only useful if 
this expansion shows up repeatedly. The result is then saved in the library in 
the installation directory, so that it can be used in future Mathematica sessions. 
Since the results are supposed to be used several times, the result of the expan- 
sion is simplified using Simplify, in order to get a more compact result. This, 
in turn, makes the extension of the library even more time consuming. 

The library files are copied to the installation directory of the package^ 
during the installation. Further libraries can be added as dcscriebed at [28]. 

Due to our naming conventions for the entries in the library files, the ex- 
pansions of HF's of parameters higher in absolute value than 9 are not allowed. 
This might be fixed if such HF's are required. 



5.4 Note on the expansion for negative parameters 

Let us consider 2Fi{—m + a, —b, —m — l + (3\ x) for m, b being positive integers, 
b> m and a, (3 small. Using the definition one gets 



n=0 



We are interested in the terms for n between m and b. All further terms vanish, 
since then {—b)n becomes 0. They are equal to 

(-m + a)...(l + a)a(-6)(l-b)...(m-6) 

~X ~r * * * ~r 



(-m + /? - l){-m + /3 + 1-0 •••(! + /? - 0(/3 - 0' 



(-m + a)...a{l + a)...{m-b + a){-b) ...{-b + m + l) 
(-m + /3-Z)...(l+/3)/? 

(105) 

If one wants to define a value for 2Fi{—m, —6, — m — l;x) one has to take the 
limit of the above expression for a and /3 going to zero. The result depends on 



^which is stored in the Mathematica variable SHypExpPath 
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the way one approches (0,0) with a and /3. In [17], one can find the formula 



2F,{-m, b, -m-l;x) = Y. /-^^^l^x" (106) 
^ (-/ - m)„n! 

whieh is also the result Mathematica gives. This corresponds to a trajectory 
in the (a, /?)-plane going along the /3 axis. Taking a trajectory along the a 
axis would lead to a 1//3 pole. Any other trajectory gives a constant times a 
function. This function happens to be the second solution of the differential 
equation satisfied by 2-^1 (—to + a, —6, —m — / + /?; x), 

x{l - x) w"{x) + {Bi - (Ai +A2 + l)x) w'{x) - Ai A2 w{x) = (107) 

with 

Ai^-m + a, A2^-b, Bi = -m - 1 + (3 . (108) 

In the case of negative or vanishing c, the value at a; = of this function is also 
equal to unity. This prevents us from discriminating the two solutions by their 
values at a; = 0. 

Since the use of the Kummer identity (|18|l 

2Fi{Ai, A2; Bi;x) ^ [1 - xf'-^'-^' 2Fi{Bi - Ai, Bi - A2; B^; x) (109) 

induces a rotation in the (a, P) plane and since Mathematica always chooses the 
trajectory along the /3 axis, the result for the HF and its Kummer transform 
will not be identical in Mathematica for this particular case. 

5.5 Performances and limitations 

The limits are set by the CPU and the amount of memory available. In all prac- 
tical cases known to the authors, however, the result is given in a reasonable 
time. The following table shows the CPU time dependence for the expansion 
of some hypergeometric functions to different orders on a 3 GHz processor/ 1.5 
GB RAM machine. 





2 


3 


4 


5 


2Fi(l + e,l-e;2 + 2e,a;) 


< 1 s 


< 1 s 


< 1 s 


7 s 


2^1(1 + aie, 1 + a2e; 2 + /3ie, x) 


< 1 s 


< 1 s 


< 1 s 


6 s 


3^2(1 + 2e, 1 - e, 2 - 3e; 1 + 3e, 2 + e,x) 


< 1 s 


< 1 s 


< 1 s 


3 s 


3^2(1 + "ic, 1 + Q;2e, 2 + 0136; 
l + /3ie,2 + /32e,a:) 


< 1 s 


< 1 s 


1.5 s 


3 s 


4^3(1 + aie, 2 + a2e, 3 + ase, 4 + a^e; 
f3ie,l + P2e,l + I33£,x) 


12 s 


20 s 


50 s 


140 s 



This package was developed in Mathematica 5 . and should work on newer 
versions. 



6 Summary 

In this paper we presented the Mathematica package HypExp for expanding ar- 
bitrary hypergeometric functions to arbitrary order in a small quantity around 
integer- valued parameters. These expansions are required for example in the 
computation of multi-loop or multi-particle phase space integrals in dimension- 
ally regularized quantum field theory. 

A first application was presented in ref . [6] , namely the e-expansion of the 
two-loop quark and gluon form factors whose exact analytic expressions contain 
hypergeometric functions in addition to F-functions. 
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A Useful relations 

In this appendix we collect useful relations among logarithms, polylogarithms 
Lin, and Nielsen polylogarithms Sn,p as well as some additional integrals. The 
package updates PolyLog as described in section 15.21 in order to allow for the 
implementation of those relations that relate polylogarithms and Nielsen poly- 
logarithms, i.e. Eq. I|113|l - l|136|l . The relations are based on [19,26] and hold 
at least for all z G W, where W is defined in Eq. ©. 



A.l Relations between logarithms and polylogarithms 



ln( 
ln( 
ln( 



z 


z 




1 




1 




1 




z 




z 




1 




z 



) = In(-z) - ln(l - z) 

) = -ln(l-z) 

) = -ln(l-^)+ln(z) 



(110) 
(111) 
(112) 



Li2(l - z 

Lt2i 



z-1 

1 



Ll2{- 



1- z 
z 

,z-l 



= -Lt2{z) + — -Hz)lnil-z) (113) 



= -Lt2iz)-^ln^{l-z) (114) 

1 TT^ 

= Li2(2)-- ln2(l-z) + — +ln(l-z) In(-z) (115) 

= - — ~ Li2(z) 

2 ^ z' 6 ' 

- xln2(--) + _-ln(-)ln(^)+Lz2(z) (117) 



(116) 



Liz{ -) = -i:7;3(z) - Lz3(l - ^) + C(3) + — ln(l - z) 

z — 1 6 

-i ln(z) ln2(l-z) + i log3(l-z) (118) 
iz3(^) = i ln3(l-z)-l In(-z) ln2(l-z) + i ln(z) ln2(l-z) 



L.3(-) 



-— - z) + Lh{l - z) 



-K- , , I. 1 , 3, 1. 



L.3(.) - y ln(--) - - ln^(--) 



(119) 
(120) 



-Li3(z) - Li3(l -z) + C(3) - \ In^(^) - ^ In(^) 

6 z 6 z 



1 



1 



+ ^ \n\l - z) + Y ln(l -z)-^ ln'(l - z) \n{z) (121) 
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= -^^\n\l-z) + \H-z)\n\l-z)-\Hz)\n\l-z) 

2 4 

+^ ln2(l - ^) + ^ - _ z) (122) 

o 45 

ii4(^) = -i^4(— ^) - :^ ln4(^) - -T In'(^) - (123) 



z ' 24 ' z ' 12 ' z ' 360 

^^5(T3^) = Y^ln^l-^)-^ln(-z)ln^(l-^) + ^ln(^)ln^(l-^) 



~ 18 ^''"^ ~ ~ 45 ^''^^ ~ + -^'^^^ ~ 



S2a{z) = ^\n\l-z)-hn{z)lnHl-z) + ^\n'{l-z) 



z , Tr"* 



-Liz{z) ln(l - z) + C(3) ln(l - z) - ^14(1 - -z) + Lii{z) 

(125) 

52,2(1-^) = -^ln^(l-^) + ^ln(^)ln3(l-^)-iln2(^)ln2(l-^) 

2 2 

ln2(l - ^) + y Hz) ln(l - z) - ln(z) Lz3(l - z) 
+Lu{l -z)- Lh{z) - Lu{^) + C(3) In(^) - ^ 

z — 1 IZU 

(126) 

^2,2(^) = ^ln^(l-.)4ln(.)ln3(l-.) + ^ln^(l-.)-^ 

-ii3(l - Z) lll(l - ^) - iZ3(^) lll(l - Z) 

+Lu{l -z) + Lu{z) + Lu{^-) 

z — 1 

(127) 

52,2(y^) = iln(z)ln^(l-z)-iln(-z)ln3(l-z)-^ln2(l-z) 

lii(-^) In(^) lii^(l -z) + ^ ln^(--2) ln^(l - -2) - Li4{z) 
-Lu{l -z)- Lu{-^) + Li3{l - z) ln(l - z) 

2 4 

+ ^ In(-z) ln(l -z) + ^- C(3) ln(l - z) 
- ln{-z)Li3{l -z) + C(3) In(-z) 

(128) 



52,3(^) = ^ln(z)ln^(l-z)-^ln^(l-^) 
1 



18 



ln\l - z) 



+ -Lh{z) \n^{l -z)+ [Lh{1 -z)- Lh{j—^)\ ln(l - z) 
C(3) 



■In (1 — z) — Ur,[z) — L«5(l — z) — Lir,[- 



-) 



+^3,2(^)+C(5) 



(129) 
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^2,3(1- 2) = —^n{z)\n\l-z)--\l^^{z)\n\l-z) + -\n\z)\n\l-z) 

2 2 

+ ^ In(^) ln2(l - ^) - ^ In'(z) ln(l - ^) - i:i4(^) ln(l - z) 
+C(3) In(^) ln(l -z) + ^ ln(l - ^) + ^ In(^) - ^ ln'(-^) 
+ i 111^(2;) ^13(1 - z)- \n{z) Lu(l - z)+ 111(2:) Lii{z) 

+ In(^) Li4(^) - 53,2(^) + 2 C(5) - ^ C(3) 

2; — 1 D 

(130) 



S2,z{^) = -^ln^(l-^)-iln(z)ln4(l-2) + ^ln3(l-^) + 2C(5) 
Li3(l - z) \v?{l -z)-\ Lisiz) ln2(l - z) - S3,2{z) 
+ iM ln2(l _ 2) + [Lz4(l - 2) + Li4(^)] ln(l - z) 

TT* 2; 

+— ln(l - ^) - 2 Li5(l -z) + Li^iz) + Lisij^) 

(131) 

S2My^J = -^ln'(l-^) + ^M-^)ln^l-^) + ^ln2(-^)ln^(l-^) 
-1 ln(-^) ln(2) ln3(l " ^) - ln'(l " ^) " ^ H-z) 
~ In^ i-z) ]n\l -z) + ^H-z) ln'(l - z) - ^-^ln\-z) 

+ - ln2(-z) ln(z) ln2(l - z) - — \n^{-z) ln(l - z) 

+^Li3(l - ^) ln(l - z)^ - ln{-z) ln(l - z) Lis{l - z) 

4 

-Lu{l - z) ln(l - 0) - Lii{-^) ln(l - ^) - ^ ln(l - z) 

1 9 

+ - In (-z) Li3(l - z) + In(-z) Li4(l - z) + In(-z) ^14(2:) 
+ In(-z) i*4(^^) + i«5(l - 2;) - 2 Li5(z) - 2 Lhij^) 

+53,2(z) + C(5)-yC(3) 



(132) 



53,2(1 -z) = -^In5(l-z) + ^ln(z)ln*(l-z)-^ln2(z)ln3(l-z) 

2 2 

ln3(l - z) + ^ ln(z) ln2(l - z) - Lz^iz) ln(l - z) 

+C(3) ln(z) ln(l -z)-^ ln(l - z) + ^ ln(z) 

- ln{z)Li4{l - z) + Li^{z) + Li^{l - z) + Li^{ -) 

z — I 



-^ + C(5)-53,(z) 



(133) 
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S3,2(^) = -\n\l-z)--Hz)\n\l-z) + -\n\l-z) + m 

+ ^ ln2(l -z)+ [Lii{^-) - Lh{z)\ ln(l - z) + 2Lh{z) 
2 ^; — 1 



ln(l -z)- Li5{l -z) + 2 Lz5(^) - S3,2{z) 
yu 2; — 1 



(134) 



^3,2(y^) = ^ln(-z)ln^(l-z)-^ln(z)ln^(l-z) 

ln^(-^) ln^(l - -2) + ^ In(^) lii^(l - z) 

+^ W(l - z) - ^ In(-z) ln^(l - + ^ ln^(l - z) 

-Li4(l - ^) ln(l -z) + Lu{z) ln(l - z) - — C(3) 

6 

~ ln(l - ^) - ^ ln(-^) + \n{-z)Lii{l - z) + 2Li,{l - z) 

-Lh{z) - Lh{^) + 53,2(2) - C(3) H-z) ln(l - z) (135) 
^; — 1 

TT^ 2Q 

53,2(-l) = Y2^^^)~32^(^^ 

There exist also relations between harmonic polylogarithms Hmi,...,mk of dif- 
ferent arguments. These are implemented in the HPL package and described in 
Ref. [27]. 

A. 2 Additional integrals 

This subsection is devoted to some additional integrals yet unknown to 
Mathematica. They are, however, not implemented in the package. 



^ ln(l - u) ln2(l - zm) tt^ , 2,, x 1 , 4, 



du^ V ^ = -- + -ln^(l-.) + -ln^(l-.) 

- ^ In(^) ln^(l -z)-2 ln^(l - z) Liziz) - [Li2iz)f 
4 ln(l - z) Li3{l -z)-2 ln(l - z) Liaiz) + 2 Lh{1 -z) + 2 Ui{z) 

+ 2Li4(^-) + 2C(3) ln(l-^) (137) 

Z — L 



I' 



1 

ln(M) ln(l - m) ln(l - zm) tt** tt^ , o,, , 1 , 4,, , 
du — ^ — ^ = 7^ + 77; In (1 -z) + w7 In (1 - z) 



90 12 ' ' 24 



- ^ In(^) ln=^(l - 2) + y Li2{z) - ^ [Li2{z)f - ln(l - z) Lh{z) 

- Lii{l -z)- Lu{z) + Liii-^) + C(3) ln(l - z) (138) 



' 2, ^ , 2 



/' , In (u) In (1 — ZU) r T- . / m9 „ , ,^ N ■■- . / ^ 

du ^^-j — ^ = -[Li2iz)]^ - 2 l7i{l - z) Li3{z) 
1-u 





+2Lz4(z) + 2C(3) In(l-z) (139) 
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/' , (u) \n(l — u) \n(l — zu) , ttt^-it / r„ „i n , •. 

du — ^-^ — !^ — !^ '- = -4HPL({3,2},^) - 10S3,2{z) 



- y Li3{z) + 2 Li2{z) Lhiz) + 6 Lir,iz) - 2 C(3) Li2{z) (140) 
^ Li2{u) _ 7r2 ln(l - z) _|_ lii(^) ln^(l ^ ^) ^ 1^(1 - z) Li2{z) 



uz 3 z 2 z z 





^ Lisjl-z) Lisjz) C(3) ^^^^^ 



^ ln(l - M^) ii2(u) TT^ TT^ ln^(l - 0) ln''(l - z) 



uz 90 ^ 3 z 2Az 



^ 5 ln(^;) ln^(l - z) ^ lii(l - zf Li2{z) _^ [Li2{z 



+ 2 



62 z 2 z 

\tl{1 - z) Uj.{l - z) - z) Liz{z) Lu{l - z) 



^--iz4(^)-C(3)^^^^ (142) 
z z z — 1 z 







ln^(«;) ln^(l - w) _ 2tt 

+ — ln(l - z) \n^{z) + - lii3(l - z) \n^{z) - \n^{l - z) In^(z) 
3 3 

- 2 ln(l - z) \t?{z) Li2{z) - 2 \n{z) Li:i{l - z) + 4 ln(l - z) ln(z) Li:i{z) 

+ 4 ln(z) Lii{l - 2) - 4 ln(z) Lii{z) - 4 ln(z) L^4(^-) + 45-3 2(2) 

2; — 1 

- iln(2)ln*(l - z) - 4C(3) ln(l - z) In(^) + 2C(3) \n^{z) (143) 
6 

^ = ln^(l -z)-- -z)- — Hz) 



+ ^ In(^) ln2(l -z) + \ 111(2) ln*(l - z) -\ H{1 - z) H{z) 

- 2 ln(2) ln^(l - 2) Li2{z) - 4 111(2) lii(l - z) Lisil -z) + 2 lii^(l - z) Lisiz) 

+ 4 ln(l - z) Lu{l - 2) + 4 ln(2) Li4(l - 2) - 4 ln(l - z) Lh{ j^) 

- 4 iz5(l - 2) - 4 ii5(2) - 4 Lhi^-) + 4 53,2(2) - 2 C(3) 111^(1 - 2) + 4 C(5) 

2 — 1 

(144) 

Z 

f ln{w)Li2{w) TT^ 2/-, N , 1 1 4/1 X / M /I N 

/ _ ^ ^ y (-^ ~ -^i + Y2 ^ ' Y ^^^^ ~ ^' 



- \ 111(2) lll^(l -Z) + \n^{z) lll^(l - 2) + 111(2) lll(l - 2) iZ2(2) 

- i [^12(2)]^ + 2 ln(2) Li3(l - 2) - 2 ln(l - 2) ^43(2) - 2 ^14(1 - z) 

+ 2 Lu{z) + 2 Lu{^-) + 2 C(3) ln(l - 2) - 2 C(3) ln(2) + ^ (145) 
2 — 1 45 
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f 27r2 

jdw [Li2{'w)Y = 6z + 6 In(l-z)-— ln(l - z) - 6 z ln(l - z) - 2 111^(1 -z) 



+ 2z ln^(l - z) + 2 ln(z) ln^(l - z) - 2 z Li2{z) + 2 ln(l - z) Li2{z) 

+ 2z ln(l - z) Li2{z) + z [Li2(z)]2 + 4Li3(l - z) - 4C(3) (146) 

f^^[lMw)]_ ^ 2HPL({3,2},z) + 4 53 2W (147) 
J w 





, \Li2(w)]^ 27r2 , , , 2 1n^(l-z) 

^ \ = — ln(l - z) + 2 ln^(l - z ^ '-+ALi2(z) 

li)^ 3 z 

+ 2 ln2(l - z) ln(z) + 2 ln(l - z) £^2(2) + 2Ml^^)^^2(^) 

z 

- [^'^^^^J' +4Li3(l-^) + 2Lz3(^)~4C(3) (148) 



/du; ,^'''^'"1 - ^ ln(l - z) - ln(z) ln2(l - z) - ln(l - z) Lz2(z) 
J (l-w)^ 3 


-2Lz3(l-z) + ii:^ + 2C(3) (149) 



^ l[i,^(,)]2+i,(i_,)i,3(,) + i^ (150) 
(1 — t(j)^ 2 1 — z 







z 

jdwU^{^^) = ^ ln2(l - ^) + ^ ln^(l -^)-\ In(^) ln'(l - ^) 


- ^ In (1 - z)2 Lz2(z) - \ [Li2{z)f - ln(l - z) £^3(1 - z) 

- ln(l - z) Li3(z) + zLi4(^^) + C(3) ln(l - z) (151) 
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